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Abstract 

We  studied  the  contribution  of  internal  energy  change  to  cathode  voltage  assuming  that  the  cathode  material  was  completely  ionic  and 
only  the  Coulomb  potential  was  effective.  From  the  internal  energy  change  caused  by  one  Li+  and  one  electron  insertion,  we  derived  a 
cathode  voltage  equation  related  to  the  Coulomb  potential,  which  we  calculated  using  a  personal  computer.  We  calculated  the  cathode 
voltage  of  an  ideal  structure  Li  YM02  (R3m)  (Li-O:  2A,  Ni-O:  2A),  which  is  a  model  compound  of  LiYNi02.  From  this  calculation,  we 
found  that  Li  VM02  has  a  two-phase  co-existence  for  0  <  x  <  0.25,  and  one  phase  for  0.25  <  x  <  1.  ©  2000  Elsevier  Science  S.A.  All 
rights  reserved. 
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1.  Introduction 

Lithium  ion  cells  [1]  are  widely  used  for  portable 
equipment  because  of  their  high  energy  density.  These 
cells  employ  lithium  intercalation  materials  as  their  anode 
and  cathode  active  materials. 

Recently,  there  have  been  several  studies  predicting  the 
cell  voltage  of  cathode  active  materials  in  lithium  ion  cells 
based  on  the  first  principles  quantum  chemistry  [2-4]. 
Deiss  et  al.  [2]  predicted  the  cell  voltage  of  LiC6/Ni02. 
They  calculated  the  average  voltage  between  the  fully 
charged  and  the  discharged  states,  neglecting  the  effect  of 
the  entropy  term. 

In  this  paper,  we  study  the  Coulomb  potential  created 
by  the  atoms  of  a  cathode  active  material  for  Li  ion  cells 
and  consider  the  possibility  of  predicting  the  voltage  from 
the  Coulomb  potential.  The  cell  voltage  is  described  ther¬ 
modynamically.  We  considered  the  contribution  of  the 
internal  energy  change  to  the  cell  voltage  ( E  vs.  Li/Li+) 
assuming  that  the  cathode  material  was  completely  ionic 
and  only  the  Coulomb  potential  was  effective  in  terms  of 
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changing  the  internal  energy.  The  process  of  evaluating  the 
internal  energy  change  involves  the  insertion  of  one  Li+ 
ion  and  one  electron,  which  are  charged  particles,  into  the 
host  matrix.  Therefore,  the  Coulomb  potential  may  have  a 
significant  effect  on  the  internal  energy  change. 


2.  Coulomb  potential 


We  assume  that  a  cathode  active  material  is  an  ionic 
crystal  and  that  the  charge  is  at  the  center  of  each  atom. 
We  ignored  any  forces  between  atoms  other  than  the 
Coulomb  force. 

Let  r, .  (A)  be  the  distance  from  a  certain  atom  k  to  the 
surrounding  atoms  of  species  i  in  an  ionic  crystal,  and  Zf 
be  the  charge  number  of  species  i.  Then  the  Coulomb 
potential  at  atom  k  is  expressed  by: 


®k  = 


e 2  Z,  Z, 

- - E- =  14.399  L- 

4nso  i,j  u  r,'j 


(1) 


where  <Pk  (eV)  is  the  Coulomb  potential  at  atom  k,  e  is 
the  elementary  charge,  s0  is  the  permittivity  of  free  space. 
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3.  The  voltage  measured  from  Li  /Li 


The  basic  equations  for  calculating  the  voltage  are 


derived  from  thermodynamics. 

A  Gr 

E  = 

nF 

(2) 

(  BG  \ 

ag'M^L 

(3) 

G  =  U  +  PV-TS 

(4) 

where  E  is  cell  voltage  (V),  A Gr  is  the  reaction  Gibbs 
function  (J /mol),  £  is  the  reaction  coordinate  (mol),  n  is 
electrons  per  molecule  oxidized  or  reduced,  F  is  the 
Faraday  constant  (C/equiv.),  G  is  the  Gibbs  free  energy 
(J),  P  is  the  pressure  (Pa),  T  is  the  temperature  (K),  U  is 
the  internal  energy  (J),  V  is  the  volume  (m3),  and  S  is  the 
entropy  (J /K). 

If  we  consider  Li  insertion  into  host  matrix  (HOST), 
AGr  is  obtained  from  the  difference  of  G  for  the  following 
reaction  divided  by  d£. 


3£Li  +  Li /HOST)  ->  Lix+a^  (HOST). 
From  Eqs.  (2)-(4),  we  have: 


1 

E= - 

nF 


dU\  (bV\  (  dS\ 

—  +  n  —  -  T  — 

J  P  T  \  d€  J  FT  \  /  FT 


(5) 


(6) 


We  will  ignore  the  second  term  in  Eq.  (6)  because  the 
value  is  very  small.  The  third  term  (Eent)  is  the  entropy 
term,  which  is  caused  by  the  number  of  different  states  of 
the  lithium  arrangement  in  the  host  matrix.  Let  N  be  the 
total  number  of  Li  sites  in  the  host  matrix,  ‘ a ’  is  the 
number  of  Li  atoms  inserted  in  the  host  matrix,  and  NA  is 
Avogadro’s  number.  The  number  of  ways  that  indistin¬ 
guishable  Li+  atoms  can  be  placed  in  N  numbered  sites, 
no  more  than  one  to  a  site,  is  N\/[a\( N  —  a)\].  If  we 
assume  an  equal  probability  of  Li  occupation,  O  = 
AM /[a\(N  —  a)\],  we  obtain: 

kT  l  Bln  LL  \ 


£“'  nF  I  3£ 


P,T 


kTNA 

nF 

-In 

RT 

=  — In 
nF 


In 


AM 


( N-a -  l)!(a+  1) 


AM 


(N  —  a)  \a\ 
N-a 


a  +  1 


RT 

nF 


In 


■  ~  A 

y 
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accurate  approximation  of  the  insertion  system  is  discussed 
by  Kudo  and  Hibino  [6].  However,  for  simplicity,  we  do 
not  adopt  their  approximation  here. 

The  first  term  in  Eq.(6)  is  the  internal  energy  term 
(E[nt).  Let  Ua  be  the  internal  energy  change  caused  by  the 
insertion  of  one  Li.  Then  we  have: 


F  = 


1  IdU 
nF  I  3£ 


P,T 


u 


ne 


(8) 


Ua  corresponds  to  the  energy  change  that  occurs  during 
the  process  in  which  Li  metal  is  separated  into  one  Li+ 
and  one  electron  followed  by  the  insertion  of  one  Li+ 
atom  and  one  electron  into  a  host  matrix. 

Let  4i  be  the  ionization  potential  (eV)  of  Li  metal,  and 
<PLi  be  the  Coulomb  potential  at  a  vacant  Li+  insertion 
site  such  as  the  site  in  Ni02.  If  we  ignore  any  contribution 
besides  the  Coulomb  potential,  the  increase  in  internal 
energy  caused  by  the  insertion  of  Li+  is  <f>Li.  Next,  we 
insert  an  electron.  We  assume  that  the  electron  is  inserted 
into  the  nearest  metal  atom  of  the  inserted  Li+  site.  The 
increase  in  internal  energy  caused  by  the  insertion  of  an 
electron  is  assumed  to  be  —  <PM  +/M,  where  <f>M  is  the 
Coulomb  potential  of  the  electron  insertion  site  and  IM  is 
the  LUMO  energy  of  an  isolated  metal  atom  in  vacuum. 
When  0M  is  calculated  for  a  cathode  before  the  Li+ 
insertion,  we  have  to  add  the  contribution  of  the  inserted 
Li+.  Thus,  we  have  to  add  —  14.399 /Ru~m  t0  where 
7?U-m  is  the  distance  between  the  inserted  Li+  and  the 
electron.  Finally,  we  obtain  the  following  equation  for  U : 


TJ  —  +  A/r 

^  a  Li  1  1 M 


4i 


14.399 


^Li-M 


(9-1) 


where  (PLl  and  (PM  are  calculated  from  the  cathode  struc¬ 


ture  before  the  Li+ 
If  we  use  the 


and  electron  insertion, 
structure  after  the  Li 


and  electron 


insertion  such  as  LiNi02  for  the  0Li  and  0M  calculation, 
the  following  equation  (Eq.  9-2)  has  to  be  used. 

14.399 


Ua  ~  ^Li  “  +  7M  _  7Li  + 


^Li-M 


(9-2) 


Finally,  the  cathode  voltage  is: 


E  ~  Eml  +  Eeni 


1 


-  <£Li  _  +  7M  _  7Li  ±  - ~ 

ne  4irsnRj: 


0  Li-M 


RT  y 
In 


nF  1  —  y 


(10) 


where  y  =  a /N.  This  is  a  well-known  result  obtained  from 
thermodynamics  [5].  In  this  paper,  we  study  the  entropy 
term  contribution  within  the  framework  of  Eq.  (7).  A  more 


where  the  sign  after  4i  is  +  if  there  is  already  Li  at  the 
center  of  k,  and  —  if  there  is  no  Li  at  the  center  of  k.  In 
the  case  of  Li  ion  cells,  n  =  1 . 


118 


J.  Yamaki  el  al.  /  Journal  of  Power  Sources  90  (2000)  116-125 


4.  Calculation  of  cathode  voltage 

First,  we  ignore  the  entropy  term  (Zsen t)  in  Eq.  (10).  We 
will  add  the  term  after  calculating  E-nV 

LiCo02  or  LiNi02  is  used  as  a  cathode  active  material 
for  Li  ion  cells.  Such  LiM02  materials  have  a  rhombohe- 
dral  R3m  a-NaFe02  layered  structure.  The  structure  is  an 
ordered  rocksalt  structure  with  edge  sharing  M06  octahe- 
dra,  which  form  M02  sandwich  sheets.  The  sheets  are 
linked  by  Li06  octahedra.  Arai  et  al.  [7]  obtained  the 
nearest  lengths  of  Ni-O,  Li-O,  Ni-Ni,  and  Li-Li  which 
are  1.957,  2.126,  2.875,  and  2.875  A,  respectively,  for 
LiNi02.  In  this  paper,  we  assume  an  ideal  LiNi02  which 
has  the  same  nearest  length  of  2.000  A  for  Ni-O  and 
Li-O.  In  this  ideal  material,  oxygen  atoms  form  regular 
tetrahedrons  linked  at  their  edges  to  form  oxygen  layers, 
and  Ni06  and  Li06  form  regular  tetrahedral  structures. 
Therefore,  the  nearest  length  of  Ni-Ni  and  Li-Li  is  1.414 
X  2  =  2.828  A. 

The  Coulomb  potential  at  atom  k  (Eq.  1)  is  calculated 
using  VIP  BASIC  (Mainstay)  for  Power  Macintosh.  The 
program  is  given  in  Appendix  A. 


4.1.  Voltage  of  Li Ni02 

From  the  calculation,  we  obtain  <PLi  =  —  18.872  eV  and 
<PNi  =  —32.1207  eV.  Using  RLi_m  =  2.828  A,  we  have, 
14. 399/RLi_Ni  =  5.091  eV.  We  obtained  an  Iu  value  of 
—  7.1076  eV  (685.78  kJ/mol)  from  the  enthalpy  change 
for  Li(cr)  — >  Li+(g)  +  e. 

We  assumed  that  the  energy  level  of  a  Ni  atom  is  the 
sum  of  the  electron  energy  level  in  vacuum  and  the 
Coulomb  potential.  Using  DV-Xa  molecular  orbital 
method  (SCAT)  [8],  we  determined  that  the  HOMO  energy 
of  [Ni06]-1  cluster  was  —0.30858  eV,  where  the  nominal 
valences  of  Ni  and  O  are  3  and  —2/3,  respectively.  The 
contribution  of  the  Coulomb  potential  of  oxygen  is 


-28.798  eV.  Therefore,  /Nl  =-  0.30858  -  28.798  = 
-29.106  (eV). 

Finally,  we  obtained  the  potential  with  the  following 
equation,  using  Eqs.  (8)  and  (9-2). 

Eint  =  18.872  -  32.120  -  7.107  +  29.106  -  5.090 

=  3.468  (V) 

4.2.  Voltage  of  Ni02 

We  assume  that  the  cathode  structure  is  not  changed  by 
the  extraction  of  Li+.  From  the  calculation,  we  had  <PU  = 
—  12.2485  eV  (empty  Li  site)  and  <Pm  =  —38.7447  eV 
(Ni4+  site).  Therefore,  using  Eqs.  (8)  and  (9-1),  we  have: 

Emt  =  12.248  -  38.744  -  7.107  +  29.106  +  5.090 

=  0.593  (V). 

However,  the  voltage  is  too  low.  Let  us  consider  a 
structure  in  which  one  Li  atom  is  inserted  into  each  plane 
along  the  c-axis.  Li+  atoms  may  occupy  the  nearest 
position  of  the  inserted  Li+  into  the  next  layer,  because  the 
Ni  valence  is  reduced  to  3,  which  decreases  the  repulsion 
force  against  Li+.  Here,  we  have  to  add  —6.73329  eV  to 
0Li  and  6.73329  eV  to  <PNi.  Using  Eqs.  (8)  and  (9-2),  we 
have: 

Eint  =  12.248  +  6.733  -  38.744  +  6.733  -  7.107 

+  29.106  -5.090  =  3.879  (V) 

Therefore,  we  have  to  use  the  following  equation  (Eq. 
11)  for  the  structure  before  Li  insertion  instead  of  Eq. 
(9-1). 

14.399 

Ua  =  +  7M  -  4i  +  -^ -  13.466.  (11) 

^Li-M 

4.3.  Lithium  ordering  in  LixNi02 

Li+  ordering  was  assumed  for  0  <  x  <  0.25,  0.25  <  x  < 
0.5,  0.5  <  x  <  0.75,  and  0.75  <  x  <  1  as  shown  in  Fig.  1. 
The  calculated  results  are  shown  in  Table  1.  Unfortunately, 


O  :  Vacant  Li  site  ^  :  Occupied  Li  site 


Fig.  1.  Lithium  ordering  assumed  for  0  <  x  <  0.25,  0.25  <  x  <  0.5,  0.5  <  x  <  0.75,  and  0.75  <  x  <  1. 
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Table  1 

Calculated  result 


x  in  Li  vNi02 

Coulomb  potential  (eV) 

Cathode  potential  (V  vs.  Li/Li  +  ) 

Occupied  Li+  site 

Ni3+  site 

Vacant  Li+  site 

Ni4+  site 

0 

- 12.2485 

-38.7447 

3.879 

0.25 

-19.4871 

-31.5061 

4.889 

0.5 

-19.5706///, 
-18.7056  = 

-31,4226///, 
-32.2876  = 

5.056///,  3.326  =  , 
average  4.191 

0.75 

-19.6541  =  , 
-18.8721/// 

-31.3391  =  , 
-32.1211/// 

5.223  =  ,  3.659///, 
average  4.441 

LOO 

-18.872 

-32.1207 

3.659 

0.25 

-11.4669  =  , 

- 12.3320/// 

-39.5263  =  , 
-38.6612/// 

2.3 14  =  ,4.044///, 
average  3.179 

0.5 

-11.5504///, 
-12.4154  = 

-39.4428///, 
-38.5778  = 

2.481///,  4.211  =  , 
average  3.346 

0.75 

-11.6339 

-39.3593 

2.648 

the  calculated  Coulomb  potential  for  two  identical  posi¬ 
tions  (shown  as  =  and  ///in  Fig.  2  and  Table  1)  did 


not  converge  to  the  same  value.  Therefore,  we  used  the 
average  value. 


x=0.5  (position  of  =) 

Fig.  2.  In-plane  position  of  Li  or  Ni,  indicated  by  “  =  ’’and  “///”  in  Tables  1-4. 
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The  equations  used  are  (for  the  structure  after  Li  inser¬ 
tion): 

E]nt  =  ~(&u-  &m)/e-  7.107  +  29.106  -5.090 

=  (-4>u  +  *m)/e  +  16-909  (12) 

where  <PLi  is  a  vacant  Li  site  and  <PNi  is  a  Ni  4+  site,  and 
(for  the  structure  before  Li  insertion): 

Eint  =  -(<£Li-  0m)/e-  7.107  +  29.106  -5.090 

+  13.466  =  (  —  <PLi  +  (PNi) /e  +  30.375  (13) 

where  <f>Li  is  an  occupied  Li  site  and  <PNi  is  a  Ni3+  site. 

There  are  four  areas  of  ordering,  which  are  0  <  x  <  0.25, 
0.25  <  x  <  0.5,  0.5  <  x  <  0.75  and  0.75  <  x  <  1.  In  each 
area,  Li+  ions  are  assumed  to  distribute  at  random.  This 
means  that  all  the  sites  are  occupied  by  Liv+  and  ey~, 
where  y  is  the  site  occupancy  in  an  area  of  ordering.  Then, 
we  have  Eint  =  E°  +  Ky,  where  E°  is  the  voltage  at  y  =  0, 
and  K  is  a  constant  (see  Appendix  B,  Eqs.  A-6  and  A-7). 
For  each  area,  we  have  to  add  the  entropy  term  Eent  = 

—  (RT/F)\n[y/(l  —  y)].  Then,  we  have: 

RT  y 

E  =  E°  +  Ky - In— —  (14) 

F  1  —  y  V  J 

where  E°  =  3.879  (V),  K  =  1.011  (V)  and  y  =  x/0.25  for 
0  <  x  <  0.25,  E°  =  3.179  (V),  K=  1.095  (V)  and  y  =  (x 

—  0.25)/0.25  for  0.25  <x<  0.50,  E°  =  3.346  (V),  K  = 
1.095  (V)  and  y  =  (x  -  0.5)/0.25  for  0.50  <  x  <  0.75, 
and  E°  =  2.648  (V),  AT  =  1.011  (V)  and  y  =  (x- 
0.75)/0.25  for  0.75  <x<  1. 

Eq.  (14)  is  shown  in  Fig.  3.  The  parameter  K  is 
positive,  which  means  that  there  is  an  attractive  interaction 
between  the  Li+  ions. 

Lithium  intercalation  compounds  are  sometimes  unsta¬ 
ble  when  their  voltage  exhibits  a  maximum  value  as  a 
function  of  the  amount  of  inserted  Li+,  as  pointed  out  by 
Armand  [5].  They  sometimes  decompose  into  two  phases 


x  in  LixNi02 

Fig.  3.  Voltage  profile  from  Eq.  (14).  (The  value  was  calculated  exactly 
using  “Mathematica”.) 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 

x  in  LixNi02 


Fig.  4.  Voltage  profile  of  LiNi02  (an  ideal  structure,  Li  ordering  was 
assumed).  The  hatched  region  is  unstable. 


similar  to  spinodal  decomposition,  even  with  an  open 
circuit,  following  the  equation: 

Li  v< HOST)  ->  !-L-EjLidy<HOST> 

/  y  —  dy  \ 

+  Li<HOST>-  (15) 

The  cathode  voltage  shown  in  Fig.  4  is  unstable  in  the 
hatched  region  and  so  may  follow  the  bold  curve  (see 
Appendix  C).  This  means  that  Li-rich  and  Li-poor  phases 
coexist.  The  plateau  voltage  is  about  0.1  V  higher  than  E° 
because  of  the  effect  of  the  entropy  term.  We  can  see  the 
phase  separation  for  0  <  x  <  0.25,  0.25  <  x  <  0.75  and 
0.75  <x<  1. 

Although  we  assumed  that  there  was  ordering,  the 
condition  for  its  appearance  have  to  be  considered.  For 
Li+  ordering  to  occur,  the  ordered  state  must  be  stable 
compared  with  the  disordered  state,  which  means  that  the 
voltage  must  be  higher  in  the  ordered  state.  We  compared 
the  voltage  as  shown  in  Fig.  5.  The  final  voltage  curve  is 
the  bold  curve  in  Fig.  5.  LiNi02  exhibits  a  two  phase 
region  for  0  <  x  <  0.25  and  is  single  phase  for  0.25  <  x  < 
1. 

Ohzuku  et  al.  [9]  reported  experimental  results  on 
LiNi02  phase  separation.  From  their  report,  the  two  phase 
region  is  0.05  <  x  <  0.25  and  other  region  is  single  phase. 
Later,  Arai  et  al.  [7]  obtained  similar  results,  and  the 
two-phase  region  is  0.15  <  x  <  0.35.  Different  results  were 
reported  by  Li  et  al.  [10].  From  their  report,  the  two  phase 
regions  are  0  <  x  <  0.3,  0.43  <  x  <  0.5  and  0.75  <  x  < 
0.85.  The  phase  separation  of  LiCo02  was  reported  by 
Ohzuku  and  Ueda  [11],  where  the  two-phase  regions  are 
0  <  x  <  0.25  and  0.75  <  x  <  0.90. 

All  the  experimental  results  show  that  the  two  phase 
region  is  0  <  x  <  0.25.  The  results  are  in  good  agreement 
with  our  calculation.  The  unstable  phases  at  x  =  0.75  and 
x  =  0.25  may  affect  the  phase  separation.  We  believe  that 
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Two  phase 


Unstable  Unstable 

two  phase  two  phase 


x  in  LixNiC>2 


unstable  unstable 

two  phase  two  phase 


Fig.  5.  Voltage  profile  of  LiNi02  (an  ideal  structure,  the  final  result). 

the  phase  separation  other  than  0  <  x  <  0.25  may  appear 
from  some  kinetic  situation  of  the  unstable  phases. 

It  is  interesting  to  consider  the  deformation  of  the 
cathode  structure  during  Li+  extraction.  The  Li+  ordering 
at  x  =  0.5  is  unstable,  however  as  mentioned  above,  it 
may  deform  the  structure  from  hexagonal  to  monoclinic. 
The  regular  hexagon  shown  in  Fig.  2  may  deform  in  one 
direction  when  ordered  Li+  atoms  are  inserted.  Thus,  the 
structure  must  be  monoclinic,  although  we  cannot  deny  the 
contribution  of  the  Jahn-Teller  effect. 

5.  Fitting  to  real  LiNi02  structure 

The  calculated  cathode  voltage  is  in  good  agreement 
with  the  experimental  value.  However,  the  structure  we 


Fig.  6.  Voltage  profile  of  LiNi02  (4%  deformed  structure  to  fit  the  real 
LiNi02). 


used  is  different  from  the  real  structure,  and  so  we  tried  to 
fit  our  structure  to  the  real  one.  As  we  mentioned,  the 
nearest  lengths  of  Ni-O  and  Li-O  are  1.957  and  2.126  A, 
respectively,  for  LiNi02.  The  mean  length  is  2.0415  A, 
and  the  Ni-O  distance  shrinks  by  4%  and  the  Li-O 
distance  expands  by  4%.  We  calculated  such  a  case,  and 
the  results  are  shown  in  Table  2.  The  equations  we  used 
for  the  calculation  are  (for  the  structure  after  Li  insertion): 

E  =  O.98(0Li  —  0Ni)/<?  —  7.107 

+  29.106-0.98  X  5.090 

=  0.98(  -  0Li  +  &m)/e  +  17.010  (16) 

where  <PLi  is  a  vacant  Li  site  and  0Ni  is  a  Ni4+  site  with 


Table  2 

4%  deformed  LiNi02 


x  in  Li A  Ni02 

Coulomb  potential  (eV) 

Cathode  voltage 

Occupied  Li+  site 

Ni3+  site 

Vacant  Li+  site 

Ni4+  site 

(V  vs.  Li/Li +  ) 

0 

-11.3905 

-39.6177 

2.545 

0.25 

- 18.6292 

-32.379 

3.535 

0.5 

-18.7127///, 

- 17.8476  =  , 
average  —18.2801 

-32.2955///, 
-33.1606*  =  , 
average  —32.7280 

average  2.851 

0.75 

-18.7961 

-18.0142/// 

-32.2121 

-32.994/// 

average  3.096 

1.00 

-18.0413 

-32.9936 

2.356 

0.25 

- 10.609  =  , 
-11.474///, 
average  —11.0415 

-40.3992  =  , 
-39.5342///, 
average  —39.9667 

average  1.860 

0.5 

- 10.6924///, 
-11.5575  =  , 
average  —11.1249 

-40.3158///, 
-39.4507  =  , 
average  —39.8832 

2.023 

0.75 

- 10.7759 

-40.2323 

1.339 
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O  :  Vacant  Li  site  ^  and  Q  :  Occupied  Li  site  of  x=0.9375 


:  Occupied  Li  site  of  x=0.75 


Fig.  7.  In-plane  position  of  Li  at  ordering  of  *  =  0.9375-. 


a  4%  deformation  from  2  A,  and  (for  the  structure  before 
Li  insertion): 

E=  — O.98(0Li-  &m)/e~  7-107  +  29.106  -  0.98 

X  5.090  +  0.98  X  6.733  X2 
=  0.98(  -  (PLl  +  0Nl) /e  +  30.207  ( 17) 

where  &>u  is  an  occupied  Li  site  and  <PNi  is  a  Ni3+  site 
with  4%  deformation  from  2 A. 

The  final  result  of  the  cathode  voltage  dependence  on 
the  amount  of  inserted  Li  is  shown  in  Fig.  6  (bold  line). 
The  voltage  becomes  lower,  and  its  difference  from  the 
experimental  value  increases.  The  shape  in  Fig.  6  is  almost 
the  same  as  that  in  Fig.  5  when  the  latter  is  shifted  1 .4  V  in 
the  lower  voltage  direction.  This  means  that  the  phase 
separation  does  not  have  any  effect  by  this  fitting  of  the 
structure.  However,  it  has  to  be  recognized  that  even  a 
small  cathode  deformation  causes  a  great  change  in  the 
voltage.  We  believe  that  we  had  a  lower  voltage  because 
the  /Ni  value  we  used  was  small.  If  we  estimate  /Ni  from  a 
(Ni06)-9  cluster,  /Ni  =  34.091  eV,  though  the  value  is 
very  large,  and  the  voltage  increases  by  about  5  V.  We 
need  a  more  precise  estimation  of  /Ni. 


6.  Other  Li  ordering 

At  x  =  0.75  +  (3/4)  X  0.25  —  0.9375,  there  is  ordering 
as  shown  in  Fig.  7.  The  results  are  shown  in  Tables  3  and 


Table  3 

4%  deformed  LiNiCL 

x  in  Li  A.Ni02 

Coulomb  potential  (eV)  Cathode  voltage 

Occupied  Li +  site 

Ni3+  site  (T  vs.  Li/Li+) 

0.9375 

-17.2896///, 

-33.7186///,  1.674 

-18.0699  =  , 

-32.9383  =  , 

average 

average 

- 17.6797 

-33.3284 

Table  4 
Ideal  LiNi02 


vinLiYNi02  Coulomb  potential  (eV)  Cathode  voltage 

Occupied  Li +  site 

Ni3+  site  (V  vs.  Li/Li +) 

0.9375  -18.1476///, 

-32.8456///,  2.991 

-18.9279  =  , 

-32.9383  =  , 

average 

average 

-18.5377 

-32.4554 

4.  The  voltage  is  very  low.  Therefore,  such  ordering  does 
not  appear. 


7.  Tetrahedral  Li  site  of  ideal  Ni02 

From  the  calculation,  we  had  <PLi  =  —  9.930  eV  (vacant 
Li  tetrahedral  site)  and  <PNi  =  —38.7447  eV  (Ni  4+  site). 
RLi_M  is  2.3804  A  (14.399 /RU_M  =  6.048  eV).  Using  Eq. 
(11),  we  have: 

E  =  9.930  +  6.790  -  38.744  +  6.790  -  7.107  +  29.106 
-6.048  =  0.718  (V). 

Therefore,  Li+  is  not  inserted  into  the  tetrahedral  sites 
but  octahedral  sites. 


8.  Conclusion 

Cathode  voltage  can  be  calculated  from  the  Coulomb 
potential  in  a  cathode  material,  when  we  consider  Li 
ordering  and  decomposition  into  two  phases.  Our  calcula¬ 
tion  predicts  a  LixNi02  phase  separation  of  0  <  x  <  0.25 
and  this  is  in  good  agreement  with  experimental  results. 
We  believe  that  there  may  be  phase  separation  in  addition 
to  0  <  x  <  0.25  as  the  result  of  some  kinetic  situation  of 
the  unstable  phases  at  around  x  =  0.5  and  0.75.  The  effect 
of  the  cathode  deformation  on  the  potential  was  very  large, 
however,  the  effect  was  very  small  as  regards  the  overall 
shape  of  the  voltage  profile. 

We  are  now  planning  to  study  the  voltage  of  spinel 
LiMn204  and  other  cathode  active  materials  using  the 
method  described  in  this  paper. 
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Appendix  A 

The  program  we  used  to  calculate  the  Li  site  potential 
of  the  ideal  LiNi02  is  given  here.  The  summation  was 
performed  for  20  X  20  X  20  =  8000  unit  cells.  We  used 
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the  potential  at  an  oxygen  atom  to  check  the  reliability  of 
the  summation.  The  oxygen  potential  is  theoretically  dou¬ 
ble  the  Madelung  constant  of  NaCl.  The  difference  be¬ 
tween  the  calculated  and  theoretical  values  is  0.333  V. 

N3  =  10 
N1  =  12*N3 
N  =  1*  N3 
N2  =  2  *  N3 
r  =  0.5  *  SQR(2) 
ex  =  0 
qx  =  0 

FOR  k  =  —N1  TO  N1 
FOR  j  =  —  N2  TO  N2 
N4  =  abs(N  —  int(abs(j /2))) 

N5  =  N4 

IF  j  mod  2=1  or  j  mod  2  =  —  1  THEN 
N4  =  N4  —  1 
N5  =  N4  +  1 
END  IF 

FOR  i  =  —  N5  TO  N4 
s  =  0 

p  =  0 

IF  j  mod  2  =  1  THEN 
s  =  r 
END  IF 

IF  j  mod  2  =  -  1  THEN 
s  =  — r 
END  IF 

IF  k  mod  3  =  1  or  k  mod  3  =  —  2  THEN 
p  =  SQR(3)*  2  *  r/3 
END  IF 

IF  k  mod  3  =  2  or  k  mod  3  =  —  1  THEN 
p=  —  SQR(3)*  2  *  r/3 
END  IF 

yl  =  (p+j*SQR(3)*r)*l 

xl  =  (s  +  i  *  2  *  r) *  1 .00 

zl  =  k* r  *  SQR(2/3) 

rl  =  SQR(xl  A2  +  yl  A2  +  zl  A2) 

q=  1 

IF  k  mod  4=1  or  k  mod  4  =  3  or  k  mod 
4  =  -  1  or  k  mod  4  =  -  3  THEN 
q=  -2 
END  IF 

IF  k  mod  4  =  2  or  k  mod  4  =  —  2  THEN 
q  =  3 
END  IF 

IF  (i  -  N4)  *  (i  +  N4)  =  0  THEN 

q  =  q/2 

END  IF 

IF  (k  —  Nl)*(k  +  Nl)  =  0  THEN 

q  =  q/2 

END  IF 

IF  (j  -  N2)  *  (j  +  N2)  =  0  THEN 

q  =  q/2 

END  IF 


IF  (i  -  N4)*(i  +  N4)  =  0  and  j  =  0  THEN 

q  =  q/2 

END  IF 

IF  i  =  0  and  j  =  0  and  k  =  0  THEN 
m  =  0 
print  i,  j,  k 
EESE 

m  =  q*  1  /rl 
END  IF 
ex  =  ex  +  m 
qx  =  qx  +  q 
NEXT  i 
NEXT  j 

IF  k  mod  10  =  0  THEN 

print  “k  =  ”;k,  “cx  =  ”;cx; 
print  “q  =  ”;qx 
END  IF 
NEXT  k 

print  “N  =  ”;N,“N1  =  ”;N1,“N2=  ”;N2,  “Final 
c  =  ”;cx 

print  “Total  q  =  ”;qx 

f=  14.399  *cx/2 

print  ‘  ‘  potential  (V)  =  ” ;  f 


Appendix  B 

We  estimate  the  Coulomb  potential  change  caused  by 
the  insertion  of  Li+  ions  and  electrons  into  the  site  shown 
in  Fig.  1.  Because  Li+  is  assumed  to  distribute  at  random, 
the  probability  of  site  occupancy  is  y.  After  averaging  the 
position  with  time,  we  can  derive  the  structure  in  which  all 
the  sites  are  occupied  and  Li  and  Ni  charge  numbers  are  y 
and  4  —  y,  respectively.  We  separate  the  Coulomb  poten¬ 
tial  into  the  contributions  from  Li,  Ni  and  O.  Using  Eq. 
(1),  we  have: 

® u,y  =  Au  +  yAu  +  4A^i  +  (4  -  y )  Am  +  AQ  (A-l) 

®M,y  =  Bu  +  yBu  +  45Ni  +  (4  -y)  Bm  +  Bo  (A~2) 

where  ^Uy  is  the  Coulomb  potential  when  the  Li  occu¬ 
pancy  is  y  at  the  Li  site,  <f>Ni  is  that  of  the  Ni  site,  ALi, 
A lj,  ANi  Ajjj  and  A0  are  the  contributions  from  Li,  Ni 
and  O  at  the  Li  site,  respectively,  and  Bu,  B^-,  Bm,  B^ 
and  BQ  are  those  at  the  Ni  site.  indicates  the  contribu¬ 
tion  from  occupied  sites,  and  XM  is  the  contribution  from 
vacant  sites  ( X :  A,  B\  M:  Li,  Ni.). 

Let  0Li  o  and  <&u  o  be  the  potential  at  y  =  0,  then  we 
have: 

=  ^n.o  (  ALi  —  ANi)  (A-3) 

^m,y  =  ^Ni,o  +  y(  Bu  ~  Bm)  •  (A-4) 

Thus,  we  have: 

^Li  ,>  _  ^Ni.y 

=  ^Li,o  —  ^Ni.o  y(  Au  —  ANi  —  Bu  +  Bm) .  (A-5) 
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We  consider  the  Li+  ordering  along  the  c-axis.  Liy+ 
ions  are  already  ordered.  Thus,  the  remaining  Li(1-y)+ 
ions  have  to  arrange  themselves  along  the  c-axis.  The 
contribution  to  Ua  is  —  13.466  (1  —  y). 

From  (Eqs.  (8),  (11)  and  (A-5),  we  derive: 

Eim  =  E«  +  yK  (A-6) 

where: 

K=  -(ALi-Am-BLi  +  Bm+  13.466)/e.  (A-7) 

For  ideal  Li  x  Ni02,  Au  =  BNi  and  AN1  =  BLi.  In  this  case, 
K=  ~[2(AU -Am)+  13.466]/c. 

As  shown  in  Fig.  1,  the  shapes  of  the  vacant  sites,  in 
each  ordering,  are  congruent.  Therefore,  Au,  ANi,  ^Li 
and  BNi  have  the  same  value  in  each  ordering  in  Fig.  1. 


x  in  Lix<HOST> 


x  in  Lix<HOST> 

Fig.  8.  The  shape  of  the  reaction  Gibbs  function  and  the  Gibbs  function 
derived  from  the  cell  voltage  with  a  maximum  point,  (a)  Cell  voltage,  (b) 
reaction  Gibbs  function,  (c)  Gibbs  function. 


Fig.  9.  Estimation  of  AG,  to  judge  the  spontaneity  of  the  phase  separa¬ 
tion. 

This  means  that  K  has  the  same  value  theoretically  for 
each  ordering. 

Appendix  C 

We  consider  a  case  where  a  cell  voltage  has  its  maxi¬ 
mum  at  point  A,  between  points  B  and  C,  as  shown  in  Fig. 
8(a).  From  Eq.  (2),  we  can  show  the  reaction  Gibbs 
function  (AGr)  (Fig.  8(b))  and  the  Gibbs  function  (G)  (Fig. 
8(c)).  When  the  cathode  material  shows  no  phase  separa¬ 
tion,  a  spontaneous  reaction  occurs  from  point  B,  through 
point  A,  to  point  C,  because  AGr  <  0.  This  is  a  cell 
discharge,  and  it  does  not  occur  when  the  cell  is  in  open 
circuit.  We  studied  the  phase  separation  of  the  cathode  at 
point  A  into  those  of  points  Al  and  A2,  which  reaction 
may  occur  even  in  open  circuit.  The  reaction  is  written  as: 

Li v(HOST>  ^Li^CHOST)  +  n2LiA2<HOST>  (A-8) 

where  rtl  +  n2=l,  x  =  nlxl+n1x2  and  xl<x<x2. 
The  reaction  can  be  considered  as  the  following  two 
reactions  occurring  at  once. 

nx Lix<HOST>  -nx(x-x\)U^  ^Li^HOST)  (A-9) 
n2LiY(HOST)  +  n2(  x2  —  x)Li  — >  n2LiA;2(HOST). 

(A-10) 

Let  AG1  be  the  change  in  the  Gibbs  function  caused  by  the 
reaction  of  Eq.  (A-9)  and  AG2  be  that  caused  by  the 


x  in  LixNi02 

Fig.  10.  Estimation  of  AG,  calculated  by  using  Eq.  (14). 
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reaction  of  Eq.  (A-10).  If  AG1  +  AG2  <  0,  the  phase 
separation  is  spontaneous. 

From  Fig.  9  and  Eq.  (1),  we  obtain: 

AG{  =nl[(x-xl)gl  +  0.5(v-vl)g2  +  5l]  (A-ll) 

AG2  =  —  n2[(x2  —  x)  gl  +  0.5(x2  —  x)  g2  +  52] . 

(A-12) 

Then,  using  the  relations  nl+n2=  1  and  x  =  nlx  1  + 
n2  x2,  we  have: 

A  Gt  =  AG1  +  A  G0  =  /?  1 5 1  —  n2S2. 

Thus,  the  phase  separation  is  spontaneous,  if  nlS  1  < 
n2S2. 

The  profile  of  Eq.  (14)  is  shown  in  Fig.  10.  Because 
nl=(x 2  —  x)/(x2—x\)  and  n2  =  (x  —  x i)/(x2  —  xl), 
n2  is  very  large  compared  with  nl.  The  actual  values  are 
nl  =  0.175  and  n2  =  0.825  from  Fig.  10.  At  the  end  point 
of  the  phase  separation,  the  values  must  be  nl  =  0.029  and 
n2  =  0.971.  And  52  is  larger  than  51  from  Fig.  10. 
Therefore,  the  phase  separation  occurs  spontaneously. 
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